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Abstract 


I/' 

The  coverage  dependence  of  the  IR  stretching  frequencies  for  dissociative 
Type  I  adsorption  of  Hj  and  Dj  on  ZnO  powders  has  been  measured  using  trans¬ 
mission  IR  spectroscopy.  By  comparing  the  frequency  shifts  observed  when 
the  isotopic  composition  of  the  adsorbate  is  varied  at  constant  total  coverage 
with  the  shifts  observed  when  the  total  coverage  of  either  pure  component  is 
varied,  we  can  separate  the  dynamic  and  static  contributions  to  the  coverage 
induced  frequency  shifts.  The  ZnH  and  ZnD  shifts  are  due  primarily  to  elec¬ 
trodynamic  interactions.  The  observed  shifts  are  in  good  agreement  with  the 
Hammaker  model  for  dynamic  dipole-dipole  interactions,  if  adsorption  is  assumed 
to  occur  on  (2  x  2)  reconstructed  ZnO  (0001)  surface  planes.  In  contrast, 

4 

the  OH  and  OD  shifts  are  due  to  electrostatic  and  inductive  interactions. 

The  electrostatic  contribution  can  be  estimated  using  a  modification  of 
Buckingham's  treatment  of  local  environment  effects.  A  qualitative  model, 
based  on  electron  localization  effects,  is  presented  to  explain  the  observed 
inductive  contribution. 

'N 


I.  Introduction 


The  adsorption  of  ^  on  ZnO  has  received  renewed  attention  in  recent 
years  because  ZnO  is  a  primary  component  in  commercial  methanol  synthesis 
catalysts  (1) .  Klier  (2)  has  suggested  that  the  role  of  the  ZnO  component 
in  these  catalysts  is  to  provide  sites  for  the  dissociative  adsorption  of  h2, 
since  H2  is  known  to  adsorb  dissociatively  on  pure  ZnO  (3-5).  The  active 

sxtes  for  this  so-called  Type  I  adsorption  of  H2  on  pure  ZnO  consist  of  pairs 
of  Zn  and  0  surface  atoms: 

H2  +  -Zn-0-  *  -Zn-O-  (1) 

However,  little  additional  information  is  available  about  the  geometry  of 
Type  I  adsorption  sites  on  ZnO. 

One  method  for  obtaining  information  about  the  geometry  of  adsorbed 
species  is  to  measure  their  vibrational  frequency  shifts  as  a  function  of  ad¬ 
sorbate  coverage  (6) .  These  frequency  shifts  contain  information  about  both 
the  separation  and  the  relative  orientation  of  adsorbed  species.  By  varying 
both  the  total  coverage  and  the  isotopic  composition  of  the  adsorbate,  it  is 
possible  to  separate  the  dynamic  dipole  contribution  to  the  frequency  shift 
from  the  contribution  due  to  static  dipole  or  chemically— induced  interactions 
(7). 

Previous  workers  have  reported  that  the  coverage-induced  ZnH  and  OH  fre¬ 
quency  shifts  on  ZnO  were  independent  of  isotopic  composition  (8).  However, 
while  performing  experiments  for  our  previously  reported  study  of  the  energe¬ 
tics  of  Type  I  H2  adsorption  (5),  we  observed  that  the  isotopic  frequency 
shift  behavior  was  somewhat  more  complex.  Moreover,  a  quantitative  discussion 
of  the  cause  of  the  frequency  shifts  for  Type  I  adsorption  has  not  been  given. 
Therefore,  the  present  investigation  was  undertaken. 
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1&2  and  pure  D£  adsorbed  on  ZnO,  and  for  adsorbate  layers  as  a  function 

of  composition.  By  comparing  the  pure-  and  mixed-  isotope  results,  we  show 
chat  the  ZnH  shift  is  primarily  caused  by  dynamic  dipole  interactions,  while 
the  OH  shift  is  mainly  due  to  static  interactions.  Possible  geometries  for 
the  Type  I  adsorption  sices  are  discussed,  and  evidence  is  given  supporting 
the  partially  reconstructed  ZnO  (0001)  surface.  Using  this  geometry,  we  show 
that  the  observed  ZnH  shift  is  described  quite  well  by  Hammaker's  model  for 
dynamic  dipole  interactions  (10) .  We  also  discuss  possible  causes  for  the 
observed  OH  shift,  and  conclude  that  the  static  shift  contains  two  contri¬ 
butions  with  opposite  signs:  An  electrostatic  effect  that  can  be  described 

2 

by  Buckingham's  treatment  of  local  environment  effects  (11),  and  an  inductive 
effect  caused  by  electron  localization  at  neighboring  adsorption  sites,  which 

produces  a  change  in  the  reference  adsorbate  bond  polarity. 


II.  Experimental 

The  IR  sample  cell  and  ion-pumped  gas  handling  system  used  in  this  work 

*¥- 

were  described  previously  (5) .  The  ZnO  samples  (Kadox-25  ,  New  Jersey  Zinc* 
Co.)  were  prepared  as  pressed  disks,  evacuated  to  673  K,  and  cooled  in  0,,  to 


improve  IR  transmission.  Adsorbate  gases  (H2,  HD,  and  ;  Matheson)  were  ^ * 

used  as  received.  cJ?o- 

Absolute  intensity  measurements  for  all  four  species  (OH,  OD,  ZnH,  ZnD) 
as  a  function  of  coverage  were  obtained  by  recording  IR  spectra  of  known 
amounts  of  or  D2  adsorbed  on  a  sample.  The  amount  of  each  gas  adsorbed 
was  determined  by  measuring  the  area  under  its  temperature"- programmed  desorp¬ 
tion  curve  (5).  The  integrated  IR  intensities  were  evaluated  approximately 
as  the  product  of  the  height  of  the  absorbance  peak  times  the  full  width  at 
half  maximum.  We  estimate  the  accuracy  of  specific  absorbance  values  deter¬ 
mined  in  this  way  to  be  -+  10%. 
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Infrared  frequency  shifts  in  the  pure  adsorbate  experiments  were  measured 


at  300  K  by  varying  the  pressure  of  or  in  the  cell  from  zero  to  100  torr. 
Adsorbate  coverages  in  these  experiments  were  determined  from  the  IR  inten¬ 
sities.  For  the  mixed  isotope  experiments,  the  sample  was  exposed  to  H0:D- 

z  z 

mixtures  of  various  compositions  at  a  constant  total  pressure  of  88  torr, 
in  order  to  maintain  saturation  coverage  of  the  Type  I  sites. 


III.  Results 


Survey  spectra  of  pure  and  pure  adsorbed  separately  on  ZnO  are  shown 

in  Fig.  1.  The  lower  curve  shows  the  background  absorbance  of  an  evacuated 

sample,  with  features  at  3700-3400  cm  3  and  1600-1300  cm  3  due  to  residual 

OH  and  CO^  species, respectively .  The  spectra  obtained  with  75  torr  of 

either  or  are  shown  superposed  in  the  upper  curve.  Where  changes  from 

the  background  spectrum  are  observed,  the  solid  line  shows  the  features 

caused  by  H2  adsorption  and  the  dashed  line  shows  features  caused  by 

adsorption.  At  full  coverage,  adsorbed  produces  peaks  at  “  3481 

cm  *  and  u)„  „  *  1712  cm  \  while  adsorbed  D_  produces  isotopically  shifted 
znH  z 

peaks  at  u)nn  *  2579  cm  3  and  uj_  _  *  1232  cm  3, 
ou  £nJJ 

The  specific  absorbance  for  each  species,  A^,  is  evaluated  from  the 

intensity  vs.  coverage  measurements.  The  intensities  for  OH  and  ZnH  were  shown 

in  our  previous  paper,  and  are  described  by  the  algebraic  expressions: 

2 


NOH=0-33  V0-10  r0H 


NZnH-°-17  W0-05  XZnH 


(2-a) 

(2-b) 


where  is  the  amount  of  each  species  present  on  the  sample  (nmole)  and 
1^  is  the  integrated  intensity  (cm  3).  The  corresponding  expressions 
for  CD  and  ZnO  determined  in  this  work  are: 


The  specific  absorbance,  A^(cm/mole),  is  Chen  given  by  (13): 

A±  -  2.303  x  106  a.  I  /Nt  (3) 

2 

where  a  is  the  cross : sectional  area  of  the  sample  disk  (0.97  cm  ).  The 
zero-coverage  limit  of  A^  is  obtained  by  combining  the  linear  portions  of  Eqs. 
(2-a  to  2-d)  with  Eq.  (3).  The  dipole  derivative,  (3y/3r)^,  which  is  necessary 
for  the  frequency  shift  calculations  described  below,  is  given  by  (13) : 

Qy/3r)2  =  (3c2mr/NAV7:)Ai  (4) 

where  c  is  the  speed  of  light,  m^  is  the  reduced  mass  of  the  oscillator,  and 

N^v  is  Avogadro’s  number.  The  values  of  (3y/3r)^  obtained  in  this  way  are  listed 

o 

in  Table  1  (see  below).  We  note  that  the  value  for  OH  (1.27  D/A)  is  comparable 

o 

to  that  reported  by  Brown  for  CH^OH  in  CCl^  (1.47  D/A)  (14);  no  comparison  is 
available  for  ZnH. 

The  coverage  dependence  of  the  IR  stretching  frequencies  observed  in  the 


pure  component  experiments  is  shown  in  Fig.  2.  The  results  for  adsorption  are 
taken  from  our  previous  work  (5) ,  while  the  results  for  were  obtained  with 
a  new  sample  having  a  slightly  higher  surface  area.  We  see  that  io011  and  <i)QD 
decrease  as  adsorbate  concentration  increases.  The  total  shifts  are  -17  cm  ^ 


for  u;_„  and  -10  cm  for  In  contrast,  u>  „  and  u>_  increase  as  coverage 

Un  UU  £nn  CttlU 

increases.  The  total  shifts  are  also  smaller,  being  5  cm  2  for  ^ZnH  an<*  cm  ^ 


for  uqD-  These  results  are  consistent  with  the  work  of  previous  authors 


(3,  4,  8). 


In  Fig.  3,  we  show  the  frequency  shifts  observed  in  the  mixed-isotope 
experiments  as  a  function  of  gas-phase  composition.  The  total  pressure  in  the 
cell  was  kept  constant,  and  the  H:D  ratio  of  the  gas  phase  was  varied.  The 
total  pressure  was  high  enough  that  the  Type  I  sites  were  always  saturated. 
Thus,  the  total  adsorbate  concentration  is  constant,  and  only  the  isotopic 
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composition  of  the  adsorbate  layer  is  changing.  In  this  work  we  are  primarily 
concerned  with  the  maximum  frequency  shifts  observed  at  the  two  composition 
extremes,  and  show  the  shifts  at  intermediate  composition  for  completeness. 


JTherefore^We  retain  gas-phase  composition  as  the  independent  variable  in  Fig. 

3,  although  we  did  observe  evidence  for  a  slight  enrichment  of  OD  and 
ZnD  in  the  adsorbed  phase,  relative  to  the  gas  phase  H:D  ratio. 

From  Fig.  3  we  see  that  coQH  and  U3QD  are  essentially  independent  of 
isotopic  composition.  For  example,  =  3486  cm  \  regardless  of  whether 

it  is  surrounded  by  OH  or  OD  neighbors'.  This  shows  that  the  coverage-induced 
frequency  shifts  of  and  are  caused  by  static  interactions  with  neighbor¬ 
ing  adsorbates.  The  shifts  in  cj  and  uj  cannot  be  due  to  dynamic  interactions, 
because  dynamic  effects,  which  are  strongest  between  isotopically  identical 
species  and  which  should  be  negligible  between  H-  and  D-  oscillators,  would 
produce  an  isotope  effect  in  the  frequency  shifts;  no  such  effect  is  seen  in 

“oh  "00- 

In  contrast,  there  Is  an  isotope  effect  in  the  ZnD  and  ZnH  frequency 
shifts:  ^nD  *ncreases  as  concentration  of  D-species  in  the  adsorbate 
layer  increases,  and  ^ZnH  *ncreases  as  the  concentration  of  H-species  increases. 

In  fact,  the  overall  frequency  shift  and  the  linear,  behavior  with  respect  to 
composition  are  nearly  the  same  as  those  observed  in  the  pure  component  experiments 
and  appear  to  be  unaffected  by  the  presence  of  the  complementary  isotope. 

This  shows  that  the  ZnD  and  ZnH  frequency  shifts  are  due  primarily  to  dynamic  inter 
actions,  and  that  purely  static  effects  cancel  each  other  or  are  absent 
between  Zn-H  species. 
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IV.  Discussion 


Coverage-induced  frequency  shifts  in  vibrational  spectra  of  adsorbates 
can  be  induced  by  either  dynamic  or  static  interactions  between  molecules. 

These  two  types  of  shift  can  be  distinguished  experimentally,  since  dynamic 
shifts  depend  only  on  the  concentration  of  isotopically  similar  neighbors, 
while  static  shifts  are  independent  of  the  isotopic  composition  of  surround¬ 
ing  neighbors  (7) . 

Dynamic  shifts  are  generally  attributed  to  electrodynamic  interactions 
between  oscillating  dipoles.  The  theoretical  model  which  describes  dynamic 
frequency  shifts  was  developed  by  Hammaker  (10) ,  and  has  subsequently  been 
discussed  by  several  authors  (15-19) .  All  of  these  discussions  have  involved 
adsorbates  on  metal  surfaces,  so  that  the  prediction  of  frequency  shifts 
required  certain  assumptions  regarding  the  effects  of  image  dipoles.  The 
present  work  allows  a  much  more  straightforward  application  of  dynamic 
coupling  theory. 

In  contrast,  static  shifts  may  be  caused  by  either  electrostatic  or 
chemical  interactions,  or  both.  In  order  to  identify  the  magnitude  of  Che 
chemically- induced  shift,  which  is  the  quantity  of  greater  interest  to  the 
catalytic  chemist,  it  is  first  necessary  to  identify  the  electrostatic 
component  of  the  observed  shift.  To  estimate  this  quantity,  we  have  modified  Bucking¬ 
ham's  (11)  treatment  of  solvent-induced  frequency  shifts  to  describe  the  effect 
of  the  electric  field  due  to  a  neighboring  dipole. 

In  what  follows,  we  shall  first  discuss  the  most  probable  adsorption 
site  geometry.  We  then  consider  the  dynamic  and  electrostatic  frequency  shifts. 

For  both  types  of  interaction,  we  will  calculate  the  shifts  predicted  for 
two  occupied  nearest  neighbor  pair  sites  and  then  multiply  these  quantities 
bv  the  surface  Madelunc  constant  (vide  infra.)  to  obtain  the  total  predicted 
shifts,  for  comparison  with  the  observed  shifts.  The  difference  between  the 


observed  static  shift  and  the  predicted  electrostatic  shift  then  represents 
the  chemically  induced  shift.  We  will  conclude  by  offering  a  qualitative 
explanation  for  the  latter  quantity,  based  on  inductive  interactions 
between  adsorbates. 

Adsorption  site  geometry.  To  perform  the  frequency  shift  calculations,  it  is 
necessary  to  specify  the  geometry  of  neighboring  adsorbates.  Boccuzzi  et  al. 

(9)  have  suggested  that  Type  I  adsorption  occurs  on  a  polar  ZnO  (0001) 

surface.  The  outermost  layer  of  this  surface  consists  of  a  hexagonal  array 

of  oxide  ions,  with  one  fourth  of  the  ions  removed  in  order  to  eliminate 

the  intrinsic  surface  dipole  layer  that  would  be  present  on  an  ideally 

terminated  surface.  One  argument  presented  by  those  authors  in  support  of 

the  ZnO  (0001)  surface  is  that  each  missing  0  anion  exposes  a  cluster  of 

three  adjacent  Zn  cations.  Experimental  evidence  for  such  clusters  of  Zn  cations  is 

given  by  C0:H2  co-adsorption  studies,  which  indicate  that  the  Zn  cation  in  a 

Type  I  adsorption  site  has  two  neighboring  exposed  Zn  cations  (9,  26). 

Alternatively,  we  note  that  a  reconstructed  ZnO  (0001)  surface,  which 
consists  of  an  hexagonal  array  of  Zn  cations  with  one  fourth  of  the  cations  missing, 
can  also  provide  clusters  of  adjacent,  exposed  Zn  cations  (c.f.  Fig.  4). 

We  believe  there  are  two  pieces  of  experimental  evidence  which  favor  the  ZnO  (0001) 
surface.  First,  the  OH  stretching  region  of  the  ZnO  background  spectrum  is 
unchanged  by  H^  adsorption.  This  suggests  that  there  are  no  surface  OH  groups 
in  the  areas  where  H^  adsorption  occurs.'  It  is  reasonable  to  expect  that  the 
Zn-rich  ZnO  (0001)  surface  may  terminate  with  bare  Zn  cations, any 
initially  present  OH  groups  being  driven  off  as  water  during  the  400  °C 
pre-treatment.  In  contrast,  one  would  expect  at  least  part  of  the  0-rich  ZnO 
(0001)  surface  to  be  terminated  by  OH  groups.  Thus,  we  can  also 

conclude  that  the  observed  hackcround  Oil  species  are  due  to  ZnO  crystal  faces  on 

which  H.,  does  not  adsorb. 
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The  second  piece  of  evidence  is  that  the  ZnH  species  show  a  dynamic  frequency 
shift,  while  the  OH  species  do  not.  As  discussed  below,  chis  suggests  that  the  ZnH 
species  are  parallel  and  vibrate  normal  to  the  surface,  while  OH  species  are 
randomly  aligned  and/or  partially  screened  from  each  other  as  a  result  of  being 
located  in  a  sub-surface  layer.  Such  ZnH  and  OH  orientations  are  more 
likely  to  be  found  on  the  ZnO  (0001)  surface  than  on  the  ZnO  (0001)  surface. 

In  Fig.  4  we  show  a  model  of  the  (2x2)  reconstructed  ZnO  (0001)  surface  in 
which  one  quarter  of  the  Zn  cations  have  been  removed.  Possible  H^  adsorption 
sites  have  been  indicated,  subject  to  the  constraint  that  no  more  than  one  adorp- 
tion  site  exists  for  each  missing  Zn  cation.  With  this  constraint,  we  note 
that  the  average  nearest  neighbor  separation  between  Type  1  sites  will  be 
6.5  X.  Also,  we  have  drawn  the  figure  showing  the  absorbed  H  species  located 
directly  above  the  Zn  or  0  ions.  It  is  possible  that  the  H  atoms  are  laterally 
displaced  somewhat,  perhaps  closer  to  the  three-fold  co-ordinated  "hollow"  sites 
between  ions. 


Dynamic  Dipole  Interactions.  In  the  presence  of  an  isotopically 

equivalent  neighbor,  the  frequency  (expressed  as  wave  number,  cm  ) 

of  a  reference  dipole  is  shifted  upward  from  its  isolated  species  value, 

ip  ,  due  to  coupling  with  the  oscillating  electric  field  produced  by  a 
o 

neighboring  dipole.  The  new  frequency,  ,  is  given  by  (10,  15): 


Ld  12  r  \2  ,  1  /3m2  1  . 

'  ii  *  *  (u<»]  4tt‘c2  \3r-j  ^  m^  ii 


(5) 


where  <{>  is  a  geometric  factor  which  accounts  for  the  separation  and 
relative  orientation  of  the  dipoles: 

-Jt 

(E, 


i  i 


ii.) 

i 


Kii 


lUil- 


(6) 


Here  E  is  Che  electric  field  induced  by  the  neighboring  dipole  at  the 
position  of  the  reference  oscillator,  and  is  the  static  dipole  moment  vector  of 
the  reference  oscillator.  Equation  (6)  is  based  on  the  point  dipole  approxima¬ 
tion  and  the  assumption  that  the  static  and  dynamic  dipoles  have  the  same 

3 

direction.  For  parallel  nearest  neighbors  =  1/R  ,  where  R  is  the 

dipole  separation  distance. 

In  Table  1  we  list  the  predicted  dynamic  shift  for  each  species  for  a 

o 

pair  of  isotopically  similar  adsorbates  separated  by  6.5  A,  the  average  nearest 
neighbor  separation  suggested  above.  The  first  two  columns  contain  the 
parameter  values  for  (3y/3r)  and  cj0  used  in  the  calculation.  The  dynamic  shifts 
predicted  by  Eq.  (5)  are  given  in  the  third  column.  Qualitative  agreement 
with  experimental  results  is  already  apparent,  in  that  the  predicted  shifts 
per  dipole  pair  are  much  greater  for  ZnH  (0.58  cm-1)  and  ZnD  (0.25  cm-1) 
than  for  OH  (0.14  cm  1)  and  OD  (0.08  cm-1). 

The  dynamic  shifts  calculated  above  for  a  single  pair  of  interacting  adsor¬ 
bates  cannot  be  compared  directly  with  the  observed  dynamic  shifts,  because 
the  latter  reflect  interactions  with  all  neighbors.  For  adsorbed  dipoles 
oriented  normal  to  the  surface,  the  ratio  of  the  total  shift  to  the 
single  nearest  neighbor  shift  is  given  by  the  Madelung  constant  for  the 
surface  dipole  lattice  sum: 

M  =  L3  I  (1/R. . ,)3  (7) 

^  i  11  » 

where  L  is  the  nearest  neighbor  separation  and  the  summation  is  made  over 
all  neighbors  R^,.  If  we  assume  as  an  upper  limit  that  every  Zn  vacancy  in 
the  reconstructed  ZnO  (0001)  surface  creates  one  adsorption  site,  then 

the  adsorption  site  lattice  has  hexagonal  symmetry,  and  we  can  use  M  -  11.034 

(27)  to  obtain  the  tot  il  nreeicc.u  mift.  These  total,  predicted  shifts  are 
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listed  in  the  fourth  column  of  Table  1.  For  comparison,  the  observed 

shifts  are  listed  in  the  final  column.  The  latter  quantities  are  the  frequency 

shifts  observed  between  the  composition  extremes  shown  in  Fig.  3.  Recall 

that  the  results  in  Fig.  3  were  obtained  at  constant  total  coverage,  so  that  no  changes 

in  static  interactions  should  be  present. 

We  see  that  the  agreement  is  quite  acceptable  for  all  species 

except  OH,  for  which  the  predicted  shift  exceeds  the  observed  shift  by  0.6  cm 

We  believe  this  indicates  the  OH  (and  0D)  species  are  not  parallel,  but  instead 

are  tilted  away  from  the  surface  normal.  The  latter  orientation  would  result  in  a 

3 

value  of  $  <  1/R^  »  which  would  reduce  the  predicted  dynamic  frequency  shift. 

A  tilted  OH  axis  would  be  expected  from  a  valence  bond  description  of  H_  . 
adsorption  on  the  ZnO  (0001)  surface. 

Static  dipole  Interactions.  The  effect  of  changes  in  local  environment 
on  the  vibrational  frequency  of  an  anharmonic  ocsillator  has  been  treated  by 
Buckinghamj who  obtained  the  following  expression  for  solvent  induced 
static  frequency  shifts,  Auk..  (11,  12): 


Aw . . S  = 
1J 


i  hc^° ! i ' I  £  . -  * 

;  1 


32U  ] 

U !  -  3a  —ill  I 


1 


(8) 


I 


where  Be  =  [h/(8ir  rn^r^  )]  is  the  rotational  constant,  z  »[(r-re)/re]  is  the 
vibrational  coordinate,  r  is  the  bond  length  of  the  reference  oscillator. 


a  *  -( ui  x  /B  )"*  is  the  first  anharmonicity  constant  of  the  potential  energy 
e  e  e 


surface  of  the  reference  oscillator  (28),  and  is  the  interaction  energy  between 


the  reference  oscillator  and  the  surrounding  molecules. 

As  above,  we  shall  apply  Eq.  (8)  to  the  case  of  a  reference  oscillator 
interacting  with  a  single  neighboring  pair  site.  The  electrostatic 


contribution  to  the  interaction  energy  is: 


V  ’hi 


in. 


(9) 


Within  the  point  dipole  approximation,  E^  is  independent  of  the  reference 
oscillator  vibrational  coordinate,  and  the  necessary  derivatives  for  use  in 
Eq.  (8)  can  be  obtained  using  a  Taylor  series  expansion  for  y^: 

y±  ■  e  [y  +  (3y/3r)  Ar  +  ...]  (10) 

—X. 

where  e  is  the  unit  vector  parallel  to  the  adsorbate  bond  axis  and  y  is 
the  dipole  moment,  including  sign.  We  can  make  no  estimate  of  the  second  order  and 
the  higher  dipole  derivatives  for  these  adsorbates,  but  expect  them  to  be  small 
relative  to  the  first  derivative  (22).  For  example,  the  absence  of  significant 
intensity  in  the  0  -*•  2  overtone  of  the  ZnH  mode  implies  a  small  anharmonici"./, 
and  hence  very  small  values  for  the  higher  order  dipole  derivatives. 


Combining  Eqs.  (8-10),  we  obtain  our  final  expression  for  the  electrostatic 


frequency  shift,  Aw„ 


3B  r  a 
e  e 


±  ©i  <Eij  *  V 


Because  the  relative  frequency  shift.  Aw  /w0,  is  nearly  independent  of 


isotropic  identity  for  static  interactions,  we  will  only  examine  the  electrostatic 


frequency  shifts  for  w  and  w  .  The  additional  parameters  needed  for  this 

ZflH  OH 

calculation  are  r  =  1.595  X  and  a  =  -2.88  for  ZnH,  r  =  0‘.970  A  and  a  *  -2.12 
e  e 

for  OH  (23).  The  static  dipole  moment  is  estimated  to  be  y  =  -0.47  D  for 
ZnH  (1  Debye  «  10  esu-cm) ,  based  on  the  percentage  of  ionic  character 
predicted  by  the  Pauling  electronegativity  difference  (24),  and  y  *■  1.59  D  for  OH, 
based  on  the  observed  dipole  moment  of  ^0  (25).  The  absolute  values  of  (3y/3r) 
were  given  in  Table  1,  and  the  signs  are  assumed  to  be  the  same  as  for  y. 

The  frequency  shifts  predicted  by  Eq .  (11)  for  each  of  the  four 
pairings  of  reference  and  neighboring  dipoles  are  listed  in  Table  II.  For 
simplicity,  we  have  again  assumed  parallel  dipoles  and  6.5  X  for  the  nearest  neighbor 
separation,  recognizing  that  the  former  assumption  will  overestimate  the 
electrostatic  shift  for  OH  species.  The  predicted  total  shift  for  each  reference 
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species,  i,  is  obtained  by  adding  the  two  shifts  due  to  both  types  of 
neighbor  species,  j,  and  multiplying  by  the  Madelung  constant  (=11.034). 
We  do  not  include  the  shift  due  to  interactions  between  the  ZnH  and  OH 


dipoles  that  occupy  the  reference  pair  site;  since  these  two  species  are 
always  present  simultaneously  (cf.  Eq.  1) ,  their  "self-interaction"  cannot 
be  resolved  experimentally.  The  predicted  electrostatic  shifts  can  be 
compared  with  the  observed  static  shifts,  given  in  the  fourth  column  of 
Table  II.  The  lat ter  are  obtained  from  the  total  shifts  observed  in 
Fig.  2  by  subtracting  the  dynamic  shifts  found  in  Fig.  3 


i 

t 

} 


4 

l 

« 


k. 


We  note  that  the  predicted  electrostatic  shift  is  positive  for  OH  and 
negative  for  ZnH.  This  can  be  explained  physically  as  follows:  The  net 

4 

dipole  moment  at  an  occupied  H2  site  is  positive  (i.e.,U0H  +  PZnH  >  ^ ' 

Therefore,  an  occupied  neighboring  site  produces  a  negative  electric  field 

at  the  reference  site.  This  produces  a  decrease  in  the  equilibrium  OH  bond  length 

and  an  increase  in  the  ZnH  bond  length,  since  (Svj/cr)^^  0  and  ( 3y / 3r ) <  ® 

(c.f.  Eq.  9).  Since  the  anharmonicity  constant  is  negative,  the  new  equilibrium  bond 
lengths  produce  a  positive  frequency  shift  in  and  a  negative  shift  in 

Chemical  (inductive)  interactions .  The  difference  between  the  observed 
static  shift  and  the  calculated  electrostatic  contribution  is  listed  in  the 
final  column  of  Table  II  for  each  reference  species.  These  differences  are  attri¬ 
buted  to  chemical  interactions  between  adsorbates.  The  chemical  shift  is  -26 
cm  1  for  OJpy  and  +22  cm  1  for 

Boccuzzi  et  al.  (8)  have  attributed  the  observed  negative  OH  frequency 
shift  to  a  decrease  in  OH  bond  polarity,  caused  by  a  depolarization  of  the 
Zr.O  surface  brought  about  by  adsorption.  While  this  would  appear  to  be 
a  reasonable  explanation  for  the  OH  shift  alone,  we  note  that  H_,  induced 
polarization  should  also  decrease  the  ZnH  bond  polarity  and  produce  a  negative 
Z: "  shift,  however,  the  results  in  Table  II  indicate  that  a  positive  chemical 
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shift  exists  for  ZnH,  which  suggests  an  increase  in  polarity  for  the  ZnH  bond. 

Instead,  we  offer  the  following  alternative  explanation  for  chemically- 
induced  frequency  shifts  based  on  arguments  analogous  to  those  for  inductive 
effects  observed  in  vibrational  spectra  of  organic  molecules  (30) . 

As  suggested  schmatically  by  Eq.  (1),  a  vacant  site  capable  of  homolytic 
Hj  dissociation  must  be  able  to  provide  two  electrons  in  order  to  complete 
both  new  adsorbate  bonds.  Thus,  we  can  consider  adsorption  to  be  an 
electron  localizing  process  with  regard  to  electron  density  initially  present 
at  a  vacant  site.  adsorption  will  decrease  the  available  electron 
density  at  neighboring  sites  (i.e.,  a  negative  a-induction  effect  (30)).  This 
will  decrease  the  polarity. of  the  OH  bond,  in  which  0  is  the  electronegative 
partner,  and  increase  the  polarity  of  the  ZnH  bond,  in  which  Zn  is  the  electro¬ 
positive  partner.  Therefore,  the  effect  of  neighboring  adsorbates  will  be  to 
produce  a  negative  inductive  contribution  to  the  OH  frequency  shift  and  a 
positive  contribution  to  the  ZnH  shift.  From  the  results  in  Table  II,  we  can 
conclude  that  for  OH  species  the  chemical  contribution  outweighs  the  electro¬ 
static,  and  a  net  negative  static  shift  is  observed.  For  ZnH,  the  chemical 
and  electrostatic  shifts  fortuitously  cancel,  so  that  no  net  static  shift  is 
seen. 

We  emphasize  that  the  proposed  electron  localization  model  offers  a  simple 
description  of  the  chemically  induced  contribution  to  the  OH  and  ZnH  static 
frequency  shifts.  The  inductive  interaction  provides  an  equally  staightforward 
explanation  for  H^-CO  frequency  interactions  in  that  co-adsorption  system,  as 
we  will  describe  in  a  subsequent  publication  (26). 

V.  Summary 

Using  pure-  and  mixed-isotope  studies  of  adsorption  of  ZnO,  we  have 
shown  that  the  coverage-induced  shifts  of  ZnH  and  Zn D  are  dynamic  in  origin, 
while  the  shifts  of  OH  and  OD  are  due  to  static  interactions.  The  dynamic 


wZnH  an<*  uZnD  s^fts  are  Ascribed  well  by  Hammaker's  treatment  of  dynamic  . 

dipole  Interactions;  a  comparison  of  the  predicted  shift  for  a  pair  of  nearest 
neighbor  adsorbates  with  the  observed  total  shift  yields  a  value  for  the  surface 
dipole  Madelung  sum  which  is  consistent  with  an  adsorption  geometry  based 
on  a  reconstructed  polar  ZnO  surface,  assuming  an  average  site  separation  distance 
of  6.5  X.  Dynamic  and  shifts  are  completely  absent.  Also  the 
absence  of  background  OH  shifts  suggests  that  adsorption  occurs  on  ZnO 
(0001)  surfaces  that  are  free  of  residual  OH  groups. 

The  observed  static  OH  shift  and  the  negligible  static  ZnH  shift 
appear  to  be  the  result  of  completing  electrostatic  and  inductive  interactions. 

For  OH,  electrostatic  interactions  produce  a  positive  contribution  while 

* 

chemical  interactions  make  a  larger,  negative  contribution.  For  ZnH,  electro¬ 
static  interactions  produce  a  negative  contribution  while  chemical 
interactions  produce  an  equal,  positive  contribution.  The  electrostatic  shifts 
can  be  predicted  using  a  modification  of  Buckingham's  treatment  of  solvent- 
induced  shifts  to  describe  the  effect  of  electric  fields  due  to  neighboring 
dipoles.  The  inductive  shifts  can  be  qualitatively  described  by  considering 
Hj  adsorption  to  be  an  electron-localizing  process  which  causes  neighboring 
OH  bond  polarities  to  decrease  and  the  ZnH  bond  polarities  to  increase  as 
coverage  increases. 
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Figure  Captions 


1.  IR  survey  spectra  of  reversible  adsorption  of  H2  and  D2  on  ZnO  at  300  K. 

Lower  spectrum:  ZnO  background  (evacuated).  Upper  spectra:  —  P„  =75  torr 
-  PQ  *  75  torr.  H2 

2.  Coverage  dependence  of  IR  stretching  frequencies  for  pure  H2  or  D2  adsorbed 
on  ZnO. 

3.  Coverage  dependence  of  IR  stretching  frequencies  for  H2-D2  mixtures  adsorbed 
on  ZnO.  T  =  300  K,  P„  +  P_  =87.5  torr. 

"2  D2 

4.  Diagram  of  (2  x  2)  reconstructed  ZnO  .(0001)  surface  showing  possible  H- 

adsorption  sites.  2 


Table  1:  Calculated  IR  frequency  shifts  based  on  dynamic  dipole  interactions 
Comparison  with  observed  shifts  for  isotopic  dilution  at  constant 
total  coverage. 

|  3p/3r  |  u)0  (^ii)*3  Total  shift  (cm 

(D/A)a  (ca~^~)  (cm  predicted*2  observed** 


Adsorbate 

Species 


Table  II:  Calculated  electrostatic  contribution  to  IR  frequency  shifts: 

comparison  with  observed  pure  component  frequency  shifts,  corrected 
for  dynamic  contributions . 


Reference 

Neighboring 

A  eS  3 

A“ij 

Total  shift  (cm 

species 

species 

(i) 

(j) 

(cm  1) 

electrostatic^ 

observed0 

inductive^ 

OH 

OH 

1.2 

ZnH 

-0.4 

Z 

(=Au,*S) 

0.8 

9.0 

-17.0  +  1.5 

i  — 

-26 

ZnH 

ZnH 

0.8 

OH 

-2.7 

z 

(=Au>*S) 

-2.0 

-22.1 

-0.0  +  1.5 

+22 

Single  neighbor  interaction,  *  6.5  A 
Given  by  (11.034)  (Aa>aS) 

Will  include  electrostatic  and  chemical  interactions 


d) 


Difference  between  observed  shift  and  calculated  electrostatic  contribution 


IR  Spectra  of  Reversible  Adsorption 
of  H2  and  D2  on  ZnO  (T  =  300  K) 


Wavenumber  (cm 


Wavenumber  (cm' 


Coverage  Dependence  of  IR  Stretching 
Frequencies  for  H2-D2  Mixtures 


Gas  Composition  (atom  %  D) 


Figure  3 


Possible  Type  I  H2  Adsorption  Pair  Site  Geometry 
on  (2x2)  Reconstructed  ZnO  (0001) 
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induced  frequency  shifts.  The  ZnH  and  ZnD  shifts  are  due  primarily  to  electro¬ 
dynamic  interactions.  The  observed  shifts  are  in  good  agreement  with  the  Hammaker 
model  for  dynamic  dipole-dipole  interactions,  if  adsorption  is  assumed  to  occur 
on  (2x2)  reconstructed  ZnO  (0001)  surface  planes.  In  contrast,  the  OH  and  00 
shifts  are  due  to  electrostatic  and  inductive  interactions.  The  electrostatic 
contribution  can  be  estimated  using  a  modification  of  Buckingham's  treatment 
of  local  environment  effects.  A  qualitative  model,  based  on  electron  locali¬ 
zation  effects  is  presented  to  explain  the  observed  inductive  contribution. 


